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harmonic traps; the prequel article treats the one, two and three particle cases. 
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1 Introduction to Part II 


This is the second in a pair of articles that classifies the symmetries of a model system of identical 
particles trapped in one-dimension and interacting via Galilean-invariant two-body interactions. The 
first article motivated the study of this system and its symmetries and it considered examples and 
applications with one, two and three particles. This article formalizes and extends these results to the 
case of N particles, with multiple examples for = 4. 

The Hamiltonian under study is 




■V 


N 


( 1 ) 


where 


N 

( 2 ) 

i=l 

is the non-interacting Hamiltonian. It is the sum of identical one-particle Hamiltonians H} that each 
include the one-particle trapping potential The interacting Hamiltonian includes the sum 

of identical pairwise interaction potentials 


N 

V^ = Y.V^J- ( 3 ) 

i<j 

The specific case of contact interactions is treated in detail, but many results hold for any Galilean 
invariant interaction potential = V‘^{\Qi — Oil)- 

The goal is to classify the symmetries of and for asymmetric, symmetric and harmonic 
traps. Two kinds of symmetries are considered: configuration space symmetries, which are realized as 
orthogonal transformations of configuration space q € and kinematic symmetries, which 

are the full group of all unitary transformations that commute with or . Despite their impor¬ 
tance and intrinsic interest, neither this article nor its prequel consider dynamical (a.k.a. spectrum¬ 
generating) symmetries. 

Key to the methods and results of this article is the assumption that each one particle Hamiltonian 
Hi has a discrete, non-degenerate energy spectrum cti = {eg, ei, e 2 ,...}. The configuration space 
symmetry Ci is either trivial Ci ~ Zi for asymmetric traps or it is parity Ci ~ 0(1) for symmetric 
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traps. The kinematic symmetry group Ki is time translation Ki ~ for asymmetric traps, parity 
and time translation Ki ~ 0(1) x Tj for symmetric traps, or the unitary group and time translation 
for the harmonic trap Ki ~ U(l) x Tj. See Sections 2 of the previous article for more details about 
one particle symmetries. 

The results for = 2 and TV = 3 in the previous article could be established by direct calculation 
and enumeration. However, that is not practical for higher N. The goal is to develop algebraic methods 
that can be implemented on a computer, but this requires a degree of formality and abstraction that 
some physicists may find unfamiliar or unappealing, i.e. the ‘Gruppenpest’. To try and bridge that 
gap, the second section of this article starts out with an overview of the permutation group and 
its representations. A particular class of representations called permutation modules turn out to be 
intuitive and useful, especially when the tools of state permutation symmetry and the double tableau 
basis are employed. This section also discussed how to apply these results to the case where the N 
identical particles are fermions or bosons and have an internal structure that does not participate in 
the interaction, like spin. 

The third section considers the non-interacting Hamiltonian and shows that the configuration 
space symmetry group and the kinematic symmetry group are always larger groups than the 
following minimal constructions 

C X and C X (4) 

In words, these symmetry groups are at least the semidirect product of the particle permutation group 
Pat acting on the direct product of N copies of the single particle symmetry group Ci or Ki. The 
previous article developed this structure by analysis of two and three particle systems, where the 
group elements can be easily explicitly enumerated. In this article, the general result is established by 
using representation theory for P^r ~ Sn applied to permutation modules. The degeneracies of the 
non-interacting spectrum must be explained by the dimensionality of the irreducible representations 
of When or is larger than the minimal construction, then there is either an emergent few- 
particle symmetry or an accidental symmetry. Table I summarizes results for 2, 3 and 4 non-interacting 


particles. 
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The fourth section shows that there is also a minimal construction for the configuration space 
symmetry group Cn and the kinematic symmetry group K^v of the interacting Hamiltonian 

Cat C Pjv X Cl and K^v C P^v x Ki. (5) 

The interactions break individual particle symmetries and the non-interacting energy levels are split 
into irreducible representations of the smaller group Kat. Table 2 summarizes these results for 2, 3 and 
4 interacting particles. This reduction of symmetry can be exploited to make approximation schemes 
like exact diagonalization scale more efficiently for the same level of accuracy and to find algebraic 
results for level-splitting under weak perturbations. 

The fifth section considers the unitary limit of the contact interaction, with Hamiltonian denoted 
In this limit, the ordering of particles becomes a good quantum number and a new kinematic 
symmetry emerges for finite-energy states 

K^^-PatxDjvxKi. (6) 

The additional symmetry Dn ~ Sjv is called ordering permutation symmetry, and it provides and 
alternate set of observables that can be used to analyze how the 7V!-fold degenerate energy levels split 
into less degenerate levels and bands in the ‘near unitarity’ limit. 

Throughout the article the question is asked, what results are universal? Specifically, when is there 
enough symmetry for the spectrum of to be calculated from the properties of the one-particle 
system? Which properties are required for which approximations, and what can be said withont specific 
knowledge of the trap shape or the interaction? How does this change with increasing particles or 
spin components? The short answer is that the non-interacting limit and the unitary limit of the 
contact interaction are algebraically universal for any N. Only the single particle spectrum cti is 
required in order for the machinery of group representation theory to construct a complete set of 
commuting observables for those two limits. For other interactions, the interplay of the trap shape, 
number of particles, and specific interaction determine whether algebraically universal expressions exist 
for properties like level splitting of multicomponent particles under weak perturbations or near-unitary 


perturbations. 
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2 The Symmetric Group 

The previous article discusses one, two and three particles and group representation theory techniques is 
employed to make well-known or intuitive results seems reasonable and inevitable. For two and three 
particles, it is relatively easy to achieve results by manual calculation or enumeration; the biggest 
composition subspaces for three particle are only six-dimensional. When we get to four particles, 
the possibility for 24-fold degenerate energy levels, even without accidental or emergent symmetries, 
encourages us to develop more sophisticated techniques. At the core of these techniques is the theory of 
the symmetric group and its representations. This section develops the necessary framework to extend 
results for two and three particles to four and more. 

The Hamiltonian is invariant under particle permutations Pw for every kind of trap and for any 
Galilean-invariant interactions. Therefore the configuration space symmetry group and the kinematic 
symmetry group have the abstract group Sw as a subgroup^ The properties of S^r and its irreducible 
representations (irreps) are well-known (c.f. [Tl[2l3|lj) and frequently applied in few-body physics. The 
first subsection establishes the local notation and definitions for S^v elements and irreps. Experts in 
the symmetric group could probably skip this subsection. Novices in group theory may find it useful to 
read subsection 1.4 from the previous paper, which introduces some notation and definitions for groups 
and their representations. The second subsection defines compositions and permutation modules, and 
the third discusses how to incorporate spin and spatial degrees of freedom by taking direct products 
of the spatial Hilbert space JC and the spin Hilbert space S. 

2.1 Sat Definitions 

Elements of Sn can be denoted by permutations p = {ii ... iw} or cycles c = (ij ■ ■ ■ k). For example, 
the permutation p = {1324} and cycle c = (23) realize the same element of S 4 . All elements with 
the same cycle structure form a conjugacy class. A partition of A^ is a set of non-negative integers 
[/r] = [piP 2 ■ ■ - Mr] that sum to N. Denote the set of partitions of N by P{N). There is a conjugacy 
^ The distinction between Pjv as a physical symmetry and Sn as the abstract group is useful because state 
permutation symmetry and ordering permutation symmetry Dn are also isomorphic to a symmetric 


group. 
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class of Sat for each partition [/x] G P{N). As an example, the five partitions and conjugacy classes for 
S 4 are: 

- [fill] = [I"*]: 4 one-cycles, i.e. the identity e = () = {1234}; 

- [211] = [21^]: 1 two-cycle and 2 one-cycles, also called transpositions, i.e. the six transpositions 
(12) = {2134}, (13), (14), (23) , (24), and (34); 

- [22] = [2^]: 2 two-cycles, i.e. the three disjoint, double transpositions (12)(34) = {2143}, (13)(24), 
and (14) (23); 

- [31]: 1 three-cycle and 1 one-cycle, i.e. the eight permutations (123) = {2314}, (132), (124), (142), 
(134), (143), (234), and (243); 

- [4]: 1 four cycle, i.e. the six permutations (1234) = {2341}, (1243), (1324), (1342), (1423), (1432). 
For the symmetric group, an element and its inverse are in the same conjugacy class. All elements in 
a conjugacy class are even or odd depending on whether they can be generated by an even or odd 
number of transpositions. 

For each partition of N there is also an irreducible representation (irrep) of Sat. These irreps form 
building blocks for other representations of Sat, like the configuration space representation 0{p) or 
the unitary Hilbert space representation U{p). Irrep labels can be depicted as Ferrers diagrams, i.e. r 
rows of boxes with the tth row having pi boxes (also called Young diagrams). The Ferrers diagrams 
for TV = 4 in order from least to greatest are 



There is a canonical ‘lexicographic’ ordering of irreps: the lowest partition [A^] is the one-dimensional, 
totally symmetric irrep and the highest partition [l'^] is the one-dimensional, totally antisymmetric 
irrep. Other partitions correspond to multi-dimensional representations with mixed symmetry under 
permutations. The notation [p]^ indicates the conjugate irrep of [p]. A conjugate irrep is the partition 
of N that has the Ferrers diagram with rows and columns reversed, e.g. [31]^ = [21^], and [iV]^ = [1^]. 
Some partitions are self-conjugate, like [2^] for TV = 4 and [31^] for TV = 5. 

The finite-dimensional vector space that carries the irrep [p] is denoted A4 . Denote the dimension 
of Sat irreps by d[p]. For TV = 4, these dimensions are (i[4] = 1, (i[31] = 3, c?[2^] = 2, d[2l‘^] = 3, and 
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^[1"^] = 1. The irrep dimension d[/i] can be calculated using the Frobenius formula, the hook-length 
method, or by counting the number of standard Young tableau>0 that are possible for a given Ferrers 
diagram, e.g. for irrep [21^] the standard Young tableaux are 


1 

2 

1 

3 


1 

3 


2 


2 

4 

1 

4 

5 

3 


The Sat irrep matrices 1^^^ depend on the basis chosen for and there are several methods for 

selecting bases and generating these matrices. This article uses a standard Young tableau Y to label a 
basis |Y) e so that p S is represented as 

U{p)\Y)=Y,d^UpW)- (7) 

Y' 

The basis can always be chose so that the matrices D_ are real and orthogonal; a standard choice is the 
Yamanouchi basis convention mm in which the subgroup chain Sjv D Sjv_i D • • • D S2 is diagonalized. 
For the totally symmetric representation [N], the ‘matrices’ = 1 are trivial for all p G Sjv- For 

the totally antisymmetric representation, the matrices have the property D^^\p) = TTp, where is 
the signature of the permutation p, i.e. whether the permutation is even TTp = 1 or odd TTp = —1. 

2.2 Compositions and Permutation Modules 

Consider a sequence of N non-negative integers n = (ni, 712, • • • , un )- The composition of n describes 
the numbers rii that appear in n and their degeneracies Vi without regard to the particular sequence. 
One notation for a composition is \ v\ = \ .. .J, dropping terms with Vi — {) and omitting 

the exponent when Vi = 1. For example, the sequence (2, 0,1,4,1) has the composition [01^24J. Instead 
of numbers, more general symbols can be used, e.g. the composition of a sequence of 5 symbols {a/3j3aj) 
is denoted . Note that Vi = N and the shape [v] of a composition [i/J must be a 

partition [v\ G P{N), e.g. if \y\ = [01^24J then [v] = [21^] and if \v\ = then [v] = [2^1]. 

The set of all sequences n with the same composition [z/J forms a basis for a representation space 
of Sat called a permutation module The action of p G Sat on a basis sequence is 

p • (ni,n2, • • • ,nN) = (upi.Up^, • • • ,np„). (8) 

^ Standard Young tableaux are filled with numbers or labels that must increase to the right and to the 


bottom. 

















The dimension d\i'\ of M^'^\ or equivalently the number of sequences with composition depends 
only on the shape [i'] of the composition The formula for is 


d[v\ = 


Nl 


iyi\L'2\ ■ ■ - Vrl 


(9) 


For example, if = [01^24J then dl^l = 60 and if \ v\ = then d\v\ = 30. 

The permutation module is built on sequences composed of a single symbol, like \y\ = . 

It is equivalent to the lowest, symmetric irrep All other permutation modules are reducible with 

respect to Sat: 


= 0 ( 10 ) 

where [^] < [v\ means that the partition [/x] is lower than or equal to the composition shape [v\ in 
lexicographic ordering and K[^][v] is the Kostka number describing the number of times the irrep [/i] 
appears in . Methods for calculating the Kostka numbers are well-established, for example, using 
characters using combinatoric methods [ 2 ], or using the intrinsic group of the composition [ 3 ]. As 
an example, for A^ = 4 there are five types of permutation modules M^’^\ with the following reductions 


mW = AdW 

^[31] ^ j^[i] 0_Ad[3il 

© Ad^^^' 

© 3Ad[^^' © 2Ad[^''l © 3Ad[^^"' © Ad^^^' (11) 


The first permutation module Afl^l is the trivial, totally symmetric representation of S 4 . The second 
module is the reducible representation of S 4 when the composition \ii\ has only one different 

symbol; it is called the defining representation of S 4 . Notice that the permutation modules and 

are not simply reducible; multiple copies of the same S 4 irrep appear. 

Different copies of the same irrep Ad [''I can be distinguished by semi-standard Weyl tableau?!^ e.g. 
for the composition \ v\ = with shape [v] = [21^], there are two Weyl tableaux W with shape 

^ Semi-standard Weyl tableaux are filled with numbers or labels that must increase to the bottom, but may 
be the same or increase to the right. 





9 


[W^] = [31] 


a 

[ml 

anrl 

a 

All 

1 


A 



( 12 ) 


corresponding to the two copies of in 1 . In this example, the two copies of that appear 
in the permutation module 1 with composition \v\ = \_a0^^\ are distinguished by how they 

transform under the exchange of the symbols a and 7 . In addition to carrying a representation of S 4 , 
the permutation module 1 carries a representation of symbol permutation symmetry (exchanging 
of a and 7 ) that is isomorphic to S 2 . 

To generalize, permutation modules carry a representation of Sjv realized by sequence per¬ 
mutations on a composition of symbols • • -J. Except when [i'] = [A^j, this represen¬ 

tation is reducible, but not necessarily simply reducible. Additionally, permutations module carries a 
representation of the symbol permutation symmetry, denoted Non-trivial symbol permutation 

symmetry occurs when there are symbols in the composition that occur the same number of times. 
For example, in the composition [i^\ = [012^4], the symbols 0, 1 and 4 all appear once, so the symbol 
permutation symmetry is 6 [;yj = S 3 . For the composition [vl = the symbol permutation 

symmetry is = S 2 because a and /3 can be exchanged. Symbol permutation symmetry can be used 
to distinguish between different copies of the same Sat irreps that appear in Denote the 

irreps of 6 [;yj by {i^}. The interplay between sequence permutation symmetry and symbol permutation 
symmetry has multiple applications in the subsequent sections. 

For any N, when the composition shape is [v] = [1^], the module 1 carries the regular repre¬ 
sentation of Sat: each irrep space appears as many times as the dimension d[ij] of the irrep and 

the total dimension of 1 is = A^!. The regular representation also carries a representation 
of symbol permutation symmetry ©|^^j that is isomorphic to Sat, i.e. {ly} = [I'^j also. This double 
symmetry plays a special role in the case of the unitary limit of the contact interaction, as shown in 
subsection 3.2.1. 

One basis for permutation modules is provided by the basis of sequences (n). Alternatively, a basis 
for the permutation module can be labeled by the set of vectors \ WY), where W € \i'\ are all the 
Weyl tableaux possible for the composition [i/J and Y e \W] are all the Young tableaux for possible for 
the Ferrers diagram with shape [W]. Elements of the sequence permutation group p G Sat are realized 
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by operators U{p) that mix basis vectors with the same W and different Y G [W]. 

U{p)\WY)= Y. D^YiPWY'). (13) 

Ye[w] 

For modules that are not simply reducible, elements of the symbol permutation group p S mix 
basis vectors with the same Y and different W G {ly}: 

U{p)\WY)= Y D\^lv'i^>W'Y). (14) 

we{iy} 

The double tableaux basis are eigenvectors of conjugacy class operators constructed from the group 
algebra of Sn and and from their canonical subgroup chains Sat D Sat-i D ••• D S 2 and 

©at 3 ©Af-i 3 D © 2 - The explicit construction of these operators and determination of their 
eigenvalues is given in 0 0 , and examples for two and three particles are in the previous article. 
Class operators built from transpositions are applied in 0li[7] to analyze the model Hamiltonian 
for multicomponent fermions and fermion-boson mixtures. They can be efficiently implemented using 
standard computational algebra programs, but their details are not required for the results of this 
article. 


2.3 Symmetrization of Identical Particles with and without Spin 

The iV-particle Hamiltonian (interacting or non-interacting) always has Pat symmetry, and so the 
spatial Hilbert space 1C can be decomposed into subspaces corresponding to irreps [p] G P{N) 

IC= 0 (15) 

Mgp(jv) 

For trapped particles, each subspace is isomorphic to an infinite tower of irreps spaces Each 

particular copy of is an energy eigenspace and there may be multiple copies of corresponding 
to the same energy. 

For one-component fermions and bosons, the spin Hilbert space S is trivial and the total Hilbert 
space is just H ^ JC. Only states in the totally symmetric subspace can be populated by one- 
component bosons, and only the totally antisymmetric subspace JC'C ^ is available for one-component 
fermions. For fermions and bosons with J > 1 components, the spin Hilbert space S ^ C'' can ‘carry’ 
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some of the symmetry or antisymmetry required for bosons or fermions and so the mixed symmetry 
subspaces of K, are relevant for identical bosons and fermions. 

One way to treat symmetrization of identical particles with internal components that do not par¬ 
ticipate in the Hamiltonian is to reduce S into Sjv irrep spaces 

5 = 0 (16) 

[ti]eP{N) 

using standard techniques (c.f. BEli) and then reduce the tensor product % = JC ® S into irreps 
using the Clebsch-Gordan series for Stv- When the internal components are spin, then J = 2s -|- 1 and 
S also carries a reducible representation of SU(2). The total spin operator = JO - Sf and total spin 
component Sz = JOi invariant under Sn and can be diagonalized along with the Sn irreps. 

Explaining this topic exceeds the ambitions of the present article, but two important results are: 

— For each irrep [v\ in the decomposition of S, there will be a single bosonic state for each copy of 

in the sector 

— For each irrep [p] in the decomposition of S there is a single fermionic state for each copy of 
in the in the sector /C[y]T. 

Examples with two and three particles were provided in the proceeding article. As an example for 
four particles, consider the case of spin-1/2 fermions. Two-component particles can only have internal 
states with at most two-row S 4 irreps, and the spin Hilbert space can be reduced in several ways: 

5 = 5^ 0 0 5^'! (17) 


where 


5 ( 4 ] 


~ i:)(«=2) 




5 ( 2 ^] 





22)(s=o) 


(18) 


where denotes SU(2) irrep spaces. This reduction can be used to find the degeneracy and spin 
possible for spatial states: for every level in there are the five states with total spin 2 ; for every 
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level in 1, there are nine spin-1 states, three which each possible z-component; and for every level 
in 1, there will be two spin-0 states. Constructing the explicit spin-spatial states in terms of the 
particle basis is a technical challenge that increases with complexity as N and J get bigger, but it is 
an algebraically solvable problem. See oiiniiii] for recent applications of these methods to trapped 
particles with spin. 

In summary, reducing the spatial Hilbert space into irreps of Stv is useful for symmetrizing identical 
particles, as well as understanding the degeneracy of energy eigenstates and how the energy levels split 
and combine as the trap and interaction are changed. The rest of this article shows how additional 
symmetries of the interaction and the trap enrich this structure. 


3 Non-Interacting Particles 

For the non-interacting iV-particle system, denote the configuration space symmetry group as C% and 
the kinematic symmetry group as K^. The total non-interacting system inherits a minimal config¬ 
uration space symmetry group and a minimal kinetic symmetry group from its construction out of 
one-particle systems: 


G%DPnX Cf ^ and K%DPnX Kf (19) 

where means the group constructed from iV-fold direct product of G with itself and the particle 
permutation group Pat acts via a semidirect product x on the abelian, normal subgroups of Cf ^ or 
by rearranging terms in the direct product. 

Before diving into representation theory, let us physically motivate this construction. In the case 
of a symmetric well, each particle’s individual parity operator Ui commutes with all the other parity 
operators and with the total Hamiltonian . Therefore the configuration space symmetry must 
at least have an abelian subgroup 0(1)’^^ generated by the N commuting, parity operators. Similarly, 
each individual particle’s Hamiltonian Hf commutes with the total Hamiltonian. This implies that 
each particle’s individual time evolution is still a good symmetry: the clocks of individual particles 
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are not synchronized unless there are interaction^ Therefore the kinematic group must at least 
have an abelian subgroup generated by each individual particle’s time evolution. Because par¬ 
ticle exchanges commute with the total Hamiltonian, Pjv must also be subgroup of both 

and K^. However particle exchanges do not commute with the single-particle parities and Hamilto¬ 


nians, e.g. (ij)ni = The semidirect product in (19) captures that structure of the minimal 

subgroups of C% and the same way the semidirect product is useful for describing groups of affine 
transformations as semidirect product of linear transformations and translations. 

In the subsections below, the groups C% and and their irreps are explored in more detail. See 
Table [2 for a summary of results for = 2, 3 and 4. 


3.1 Configuration Space Symmetry Group 

The configuration space of N particles in one dimension is isomorphic to one particle in 

N dimensions, and therefore low N situations can be visualized and described using the terms and 
techniques of familiar geometry. Additionally, for asymmetric and for most symmetric wells, the con¬ 
figuration space symmetry group is a finite-order point group. Finite-order point groups are subgroups 
of the orthogonal transformations 0{N) on . Point groups in all dimensions are completely char¬ 
acterized and classified. Point groups in two and three dimensions are familiar to many physicists 
from applications in chemical and solid state physics and Schonflies notation is standard. There are 
several notations for extensions of these groups to higher dimensions, but Coxeter notation m is most 
convenient for my purposes because of its connection to the theory of Weyl groups and Lie algebras, 
which finds application in the closely related theory of the Bethe ansatz solutions [Hill]. 

For the asymmetric well, the one-particle configuration symmetry group is trivial Ci ~ Zi and the 
configuration space symmetry is isomorphic to the permutation group Sat. Each particle permu¬ 
tation is realized by a geometrical transformation of configuration space , for example: 

^ This observation may seem obvious, but it was a remark to this effect in m discussing the similarly 
separable problem of a particle in a cubic box that led me to understand the connection between permutation 
modules and K%. 
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— Two-cycles (ij) are reflections across the (TV—l)-dimensional hyperplane C dehned by 

Qi = Qj- 

— Three-cycles (ijk) are generated by two overlapping two-cycles {ij){jk). They are realized by simple 
rotations by ±27r/3 in the plane perpendicular to the (A^ — 2 )-dimensional hyperplane Vij n Vjk- 

— Double two-cycles {ij){kl) are generated by two non-overlapping two-cycles {ij){kl). These are 
double reflections across orthogonal hyperplanes Vij and Vki, and they are equivalent to a simple 
rotation by ± 7 r in the plane perpendicular to the (TV—2)-dimensional hyperplane Vij H Vki- 

— Four-cycles (ijkl) are simple rotoreflections (or improper reflections). They are realized by a rotation 
by ± 7 r /2 in the plane perpendicular to Vik H Vji, followed by reflection across the same plane. 


Longer cycle structures correspond to elements of equivalence classes of higher-dimensional orthogonal 
transformations, such as compound reflections, compound rotations and compound rotoreflections. 
One way to derive these properties is to construct the defining representation of Sjv on the set of basis 
vectors {qi,q 2 , ■ ■ ■, (Jn}- For example, the matrix representing the element (12345) = {23451} € S 5 is 


0(12345) 


'^0 1 0 0 0 ^ 
0 0 10 0 
0 0 0 1 0 
0 0 0 0 1 
yl 0 0 0 Oy 


( 20 ) 


The eigenvalues of 0(12345) are (1, exp(± 47 rz/ 5 ), exp(±27ri/5)}, corresponding to a compound double 
rotation of 47 r /5 and 2tt/5 in two orthogonal planes. The eigenvalue 1 corresponds to the eigenvector 
q = (1,1,1,1, l)/-\/5 which remains invariant under all permutations of five particles. 

The geometrical realization of P^v is equivalent to the point symmetry group of a regular TV-simplex: 
the digon for N = 2, the triangle for iV = 3, the tetrahedron for TV = 4, the pentachoron for TV = 5, 
etc. These are the finite Coxeter reflection groups Ajv_i. Note also that conjugacy classes of Pat are 
associated to equivalence classes of 0(TV—1), e.g. elements of P 4 ~ S 4 are reflections, simple rotations, 
and simple rotoreflections, and these are the three equivalence classes of 0(3). This restriction to 
equivalence classes of 0(TV—1) is because the one-dimensional manifold defined hy qi = q 2 = ■ ■ ■ = qN 
is invariant under the geometric realization of every p gVn- 
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For a symmetric well, the configuration space symmetry must contain Sjv x 0(1)^^, the 
semidirect product of the symmetric group Sat on the iV-fold tensor product of the reflection group 
0(1). Each individual particle parity operator 77* is realized by reflections across the (A^—l)-dimensional 
hyperplane given by qi = 0. The total parity 77 is the product of all the individual parities 


77 = 77i772 • ■ • Iljsf. 


( 21 ) 


The group Sn x 0(1)^'^ has order 2^Nl and it is also known as the hyperoctahedral group [l6l ITT] . 
These are the finite Coxeter reflection groups labeled BGn- The N = 2 and N = 3 dimensional 
examples are the point groups of a square and cube, respectively, and their irreps (and the reduction 
of those irreps by the subgroup Sat) are discussed in the previous article and summarized in Table 
The higher-dimensional point groups for the hypercubes are less well-known in physics but their irreps 
can be found by induction from the normal subgroup 0(1)^^. That is how the irrep dimension and 
multiplicity is calculated for TV = 4 in Table [T] 


For a harmonic well, there is an emergent multi-particle symmetry and is larger than minimal 


symmetry inherited from the construction (19). There is full rotational and reflectional symmetry in 
and so ^ 0{N) is the full orthogonal group in TV dimensions. The hyperspherical representation 
of this group is well-known: irreps are labeled by A and have the dimension |18j 


(TV + 2A-2)(TV + A-3)! 
diO{N),X) - x\iN-2)\ 


( 22 ) 


for TV > 2. This formula gives the familiar results <7(0(3); A) = 2A -I-1 and <7(0(4); A) = (A -f 1)^. This 
case has been examined in more detail in [19] and further applications are in preparation. 

Note that if were the only symmetry of the non-interacting Hamiltonian , then we would 
expect that the energy levels would have degeneracies corresponding to the dimensions of those irreps. 
However, they certainly do not and the explanation of spectral degeneracies requires the consideration 
of the kinematic symmetry group K^. 
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3.2 Kinematic Symmetry Group 

Denote the non-interacting 7V-particle spectrum by = {Eq.Ex.E^,- ■ ■}, which is still discrete but 
no longer non-degenerate. The spatial Hilbert space is decomposable into energy eigenspaces 

/C = 0/C'= (23) 

k 

where = E^IC’^. If the kinematic symmetry group has been completely and correctly 

identified, then each energy eigenspace JC^ will carry an irreducible representation of K)(f. In that case, 
the dimensions of the irreps of correspond to the degeneracies of energies Ej^ € cr^. This section 
demonstrates that unless there are accidental degeneracies or emergent few-body symmetries, then the 
minimal kinematic group Pjv x is sufficient to explain the degeneracies of Ej^ G cr^. The irreps 
of Pat >< are isomorphic to the permutation modules described in Section 2.2 and their reduction 
into Sat irreps is algebraically solvable in terms of one-particle observables and particle permutation 
operators. 

The elements of cr^ and their degeneracies can be determined by forming compositions of one- 
particle energies e„ in the single particle spectrum cti. The energy level with composition [i>\ = 
.. .J has energy E^^^ = voCq + i^iei -!-•••. Only a partial ordering of is possible unless the 
specific one-particle energies e„ € CTi are known; see previous article for examples. 

The tensor product of N non-interacting basis states is compactly denoted by 

|n) = Im) G |n2) G • • • G |nAr), (24) 

or alternatively \a/3 • • •). Call this basis, where each particle has a definite state, the “particle basis”. 
Note that the particle basis wave functions 

N 

'?n(q) = (q|n) (25) 

i=l 

can always be chosen as real functions for the trapped system. Each N-particle tensor product basis 
vector |n) is an eigenvector of with energy E^,^^ given by the composition \ v\ of n. The degeneracy 
of the energy level is the number of particle basis vectors with that composition. As explained 
in the previous section, this degeneracy is determined by the shape [v] of the composition [i/J. 
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The spatial Hilbert space is also decomposable into subspaces 1C spanned by particle basis vectors 
with composition 

/C = 0/CL'^J. (26) 


Note that this decomposition is not the same as (16). The sectors are infinite towers of Stv irrep 
spaces each with dimension whereas are the composition subspaces with dimension 

d[^J. Each is an irrep of Sjv X and isomorphic to a permutation module This is 

demonstrated in the next subsection by explicit construction of the d \v\ -dimensional irrep of S at x Tj^ ^. 


3.2.1 Irreps ofS^tK 

An element (p, t) of x T^^^ is a pair formed by a permutation p £ and a real fV-tuple t = 
.. ■ ,tN) The group multiplication rule in Sat x Tfis 

(p';t') • (p;t) = (p'p;t + O^(p)t') (27) 

so that 

(p,t) = (p,0) • (e,t), (28) 

where e is the identity in Sat and 0 is the identity in Tf'^. The N x N matrix 0{p) permutes the 
components of an A^-dimensional vector x: 


0(p)x = 0{p) {xi,X 2 ,...,Xn) 


— \^Pl 1 ^P2 "I ' ' • 1 ^PN 


)• 


(29) 


The matrix (p) = O ^{p) = 0{p ^) is orthogonal and all matrix elements are zero except a single 
one in each row and columrQ 


The choice of the transpose matrix as the automorphism on T^^^ in the multiplication rule (27) is 


so that the unitary representation U(p,t) of Sat x T^ has a natural realization on the particle basis 
® The representation O is another example of the defining representation of Sat - By a similarity transformation 
JjOJ_ it can be decomposed into Sn irreps [N] © [(A —1) 1]. The set of similarity transformations J that 
reduce the defining representation is the equivalence class of normalized Jacobi coordinate systems for N 
one-dimensional particles. 
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U{p,0)\ni 712,..., nAr) = \np^,np^,... ,np„) = \0{p)n) 

C/(e,t)|n) = exp(—ie • t)|n) 

U{p, t)\n) = exp(-ie • t)|0(p)n), (30) 


where e = {( 


-ni 5 ^712 5 ■ • ■ 5 ^riN 


}■ 


The orbit of any N particle basis vector |n) with composition [z/J under the representation (301 
spans the composition space . The unitary matrix representation (p, t) on defined by 


C/(p,t)|n)= i2!vi(p,t)|n') (31) 

n'G(n) 

is dimensional and irreducible with respect to Sat x T^^. The action of the subgroup D}-''^ (p, 0) 
on /CM is unitarily equivalent to the action of the symmetric group on permutation module 
Therefore, irreps of Sat x Tf ^ fall into equivalence classes, one for each shape [v] G P{N). 

An alternate derivation of the irreps of Sa7 x T^^^ uses the technique of induced representations: 
the character of the normal, abelian subgroup Tf'^ is the exponential exp(—le-t) and it is determined 
by the sequence of energies e, or equivalently by the sequence n. The distinct orbits of n under 
Sat are classified by the shape of the composition [i/]. The ‘little group’ of a composition [z/J, i.e. 
the group of transformations that leaves a canonical representative n G \y\ invariant, is = 
Si/i X S;,2 X • • • X S,y,,. The degeneracy d[z/J is the order of the coset SAr/S[j,j and there is a basis vector 
of the irrep for each element of the coset. For example, when [z^J = [012^4J, choose the normal-ordered 
representative sequence n = (0,1, 2, 2,4). The partition [u] = [21^] has little group S[ 2 i 3 ] ~ S 2 with two 
elements e and (34), corresponding to exchange of the identical symbols. The construction or irreps by 
induced representations guarantees irreducibility. Note that is not the same thing as the symbol 
permutation group, which would be S[ 012 = 4 ] ~ S 3 . 


3.2.2 The Double Tableaux Basis for Eigenstates of Hq 

The previous subsection establishes that composition spaces K. are irreps of Sat x Tf and isomorphic 
to permutation modules The particle basis |n) provides a natural basis for this group that 
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diagonalizes the time translation subgroup However, for several practical reasons, the reduction 

of by the Sn subgroup is often useful: 


— This reduction allows the use of the double tableau basis \WY). The irreps that occur in are 

labeled by Weyl tableaux W and the shapes of the tableaux \W] indicate the Sn irrep The 

specific Young tableau Y £ [IT] labels the irrep basis. See Ref. [S] for a detailed derivation of this 
for any N and explicit construction of the complete set of commuting observables for this basis out 
of conjugacy class operators for canonical and non-canonical subgroup chains. The previous article 
also includes an extended example with N = S. 

— For spinless particles, the reduction of into irreps of Sat provides a method for handling 

identical particle symmetrization. For particles with spin, the reduction of the spatial Hilbert space 
into Sat sectors provides a methods for calculating degeneracies, observables, and basis vectors, c.f. 
Section 2.3. 

— When interactions are added to the Hamiltonian, the kinematic symmetry Sat X Tf^ will be 
broken, but the symmetry will remain. One consequence is that (unless there are accidental or 
emergent symmetries) the degenerate energy levels of each composition space will split into 
levels, specihcally, one (i[/r]-degenerate level for each copy of in ^ 

— Also, when interactions are incorporated there are only matrix elements between non-interacting 
basis vectors within the same Sat irrep. Therefore, exact diagonalization in the double tableau basis 
requires fewer matrix elements to achieve the same accuracy, and the energy spectra are more easily 
interpreted in terms of spectroscopy and selection rules, c.f. Section 4.4. 
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As an example of the double tableaux basis with = 4, consider the reduction of a composition 
space from each equivalence class of S 4 ix Tf ^ irreps: 


^ ^ ^IP ^ ^11 

^La"/37j = © /C® 


)/C^ 


K^laMS] ^ ^0203 ^ ^ ^ }(jS^ ^ ^ j^\ 


)/C^ ©/c^ ©/c“ ©/c 


(32) 


Each composition space /CL© is isomorphic to the permutation module and each subspace K,^ is 
isomorphic to the S 4 irrep 

Depending on the shape of the composition, the composition space 1C also carries a representation 
of another group, called the state permutation group First consider an example with 

This space is d[i^\ = 12 dimensional and has basis vectors (here in the particle basis): 


jaa/lq), \ aaj/3), ja/Saq), \ aPja), jaqa/?), jaq/la), 

l/Saaq), j/Saqa), \j3jaa), jqaa/S), jqa/la), | 7 / 3 aa). (33) 


The action of the particle permutation group S 4 on this basis is given by (30). However, there is 
another symmetry that leaves the space IC^°‘ invariant: switching /3 and 7 . Denote the operator 


that switches these states as U{P^). This operator shuffles the basis vectors (33) into 


|aa 7 / 3 ), \ aaj3^), |a 7 a/?), |a 7 /la), \ af3a'y), \ afi^a), 

\^aaj3), | 7 Q;,da), \^(3aa), \j3aa'y), \j3a'^a), \j3^aa). (34) 


The operator U{j3^) generates a group isomorphic to S 2 but distinct from any S 2 subgroup of the 
particle permutation group P 4 . Call this group Every irrep of S 4 >< Tf^ with shape [21^] 

carries a twelve-dimensional reducible representation of fp y^^. Note that the fp |^„ 2 symmetry is 
not a symmetry of the full Hamiltonian and most composition subspaces are not invariant under 
TJ{f5^). However, it is useful for analyzing specific composition spaces like Another way to 
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say this is that U{/3j) and commute when restricted to the spaces like , even though the 

operators defined on the entire spatial Hilbert space K do not. 

Referring back to Section 2.2, we see that particle permutation symmetry P^v places the role 
of sequence permutation symmetry and state permutation symmetry ^yi,\ plays the role of symbol 
permutation symmetry. The representation of Stv ix Tf ^ is irreducible on a composition space ~ 
The operators U{p) representing the particle permutation p S Pat mix double tableau basis vector 
\WY) with the same W and different V and the operators 17(p) representing the particle permutation 
p € mix basis vectors with the same Y and different W. 

The coefficients {n\WY) that connect the double tableau basis \WY) and the particle basis |n) 
can be found using a variety of methods. In the previous article, I follow the conventions and methods 
of Ref. [3] for two and three particles. These methods can be extended to more particles, and although 
it becomes more complicated as N grows larger, it is still algebraically solvable for any N. Explicit 
values for {'n\WY) are not necessary for the results that follow. 

Note that only composition spaces with [v] = [1^] are isomorphic the permutation module 
1 and contain a totally antisymmetric irrep of Pat. Therefore, only for compositions of N distinct 
states does have an energy level that can be populated by one-component fermions. Because of 
this spatial antisymmetry, compositions with [i/] = [ 1 ^] play an important role in the case of contact 
interactions, as discussed below. 

3.2.3 Bosonic and Fermionic Spectral Isomorphism 

Select the unique sequence n in the composition \ v\ with state labels arranged in increasing order. For 
example, this is the sequence n = (0,1,1, 2,4) from [i/J = [01^24J. By adding 0 to the first element of 
the sequence, 1 to the next element and so on, a new composition (y') is produced which has the shape 
[u'] = e.g. the sequence n = (0,1,1, 2,4) with shape [21^] becomes the sequence n' = (0, 2, 3, 5,8) 
which has shape [1®]. This means that for every totally symmetric state in KX''^ (because there is always 
one copy of in every composition subspace), there is a partner antisymmetric state in J, and 

the relationship is one-to-one. This mapping is a generalization of the result of Crescimanno |20j that 
the bosonic and fermionic non-interacting spectrum have the same structure for the harmonic oscillator. 
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just shifted by a constant value. It is also related to the famous Girardeau fermionization mapping of 
identical one-component bosons m, as discussed in Section 5. 

3.2.4 The Group for Symmetric Traps 


For the general asymmetric trap, the kinematic symmetry is the minimal symmetry = Pn X 
and irreps are labeled by compositions \v\. These irreps sort into equivalence classes irreps based 
on their shapes [v\. Symmetric traps have the larger one-particle kinematic symmetry group Ki = 
0 (1) X Tj. Irreps are still characterized by compositions, but parity adds an additional quantum 
number. Even though two compositions have the same shape, they may not transform the same way 
under the single-particle parities TIi. As a result, not all subspaces with the same composition 
shape [v\ transform fall into the same equivalence class. For four particles, the possible equivalence 
class are (organized by shape) 


[4]: L+"J,L-"J 

[31]: [+3+2j,L+3-J,[-3+J,L-3-2j 

[2^]: L+?+2j,L+^-"J,L-?-2j 

[ 21 ^]: i+l +2 -1-3J 5 L+i +2 ~J) L+^ “1 “2], + 1 +2J j L~i “2 +J, L~i “2 —3J 

[ 1 "^] : L+i +2 +3+4J, L +1 +2 +3—J) L +1 +2 —1—2J) L+i “1 “2—3]) L~i “2 —3—4J • 

The notation ^ -ki -I- 2 J identihes an equivalence class of compositions with three states, one with 
negative parity and two distinct positive parity states, for example. 

The group is an abelian normal subgroup of K^, and so the induced representation construc¬ 
tion described at the end of Section 3.2.1 can also be used to identify these equivalence classes for 
Sat X (0(1) X Tt)^^. The reduction of composition subspaces into S^r irreps is unchanged. Note 
that the total number of distinct irreps of the minimally-constructed groups and are the same 
although the distribution of irrep dimensions is not. 


3.2.5 Emergent and Accidental Symmetries 

When is equal to the minimal kinematic symmetry group Pjv x each composition subspace 

is uniquely related to an energy level G a^. The reduction of the spatial Hilbert space into 
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irreps of is the same as the decomposition into composition subspaces (26). However, if there are 


coincident energy levels, i.e. two or more compositions lead to the same total energy, then there are 
additional kinematic symmetries. 

The first case is accidental degeneracies, like the Pythagorean degeneracies that occur in the infinite 
square well and its higher dimensional generalizations |12l [22] . These accidental degeneracies can be 
formulated as an ad hoc kinematic symmetry by defining operators that act as the identity in most 
energy subspaces but act unitarily in the accidentally-degenerate composition subspaces. In such a 
formulation, each accidental degeneracy requires the addition of new operators that commute with 
the Hamiltonian. Thinking of this as a kinematic symmetry is therefore not productive because these 
symmetry operators must be inferred from the degeneracies and not the other way around. Accidental 
degeneracies of this sort will not be considered further here. 

The other reason for coincident energy levels is that there is an emergent few-body symmetry, i.e. 
a symmetry that cannot be generated by the one-particle symmetries and particle permutations. The 
harmonic well is the most famous example. Its energy levels have a degeneracy larger than can be 
explained by Pn ix For the energy level huj{X + N/2) with total excitation X = the 

degeneracy is [121IM] 


(X + N - 1)1 

d(U(iV),A) = L^±^^. (35) 

This degeneracy can be derived from combinatorics, or it can be explained by the fact that ^ U(fV), 
the group of unitary transformations in N dimensions. The U(A^) symmetry can be thought of either 
as the group of symplectic, orthogonal transformations in phase space or as the group of unitary 
transformations of the N annihilation operators: Sp(2iV) n 0{2N) ~ U(iV). The irreps of U(fV) can be 
labeled by total excitation X and the irrep space is the direct sum of all composition spaces with the 
same total excitation. 


4 Interacting Particles 

The introduction of Galilean-invariant two-body interactions Vij among the identical particles breaks 
the symmetry encapsulated by the subgroups Ci^ and The levels in the non-interacting energy 

spectrum therefore split and degeneracies are reduced when two-body interactions are turned on. 
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However, the permutation symmetry subgroup is preserved, as well as any transformation in or 
Kfthat also commutes with the interaction operator 

N 

= (36) 

i<j 

Denote the symmetry groups of = Hq + by Cat and Kjv. They are subgroups of and K^, 
respectively. See Table for information about C m and Kjv and their irreps for low particle numbers. 
Before classifying them for the three types of traps, we take a brief detour into the symmetries of the 
two-body matrix elements. 


4.1 Symmetries of the Two-Body Matrix Elements 


The two-body matrix elements are the matrix elements of the two-body interaction in the non¬ 
interacting particle basis: 


(n|Ei2|n') = (nin2|Vi2|n;n2)(5„3„' • ■ ■ 


rianj n^n'^ 


(37) 


Because of the hermiticity of Vij , and remembering that the stationary states of single trapped particles 
can always be chosen as real, the two-body matrix elements are also real and have the property 


0.3 - 7(5 

V = ^«/3' 


(38) 


Galilean invariance constrains the position representation of the two-body interaction to have the form 


= V ih -qj\)S (q-q')- 


(39) 


where is a scalar function of particle separation. The potential V'^ could be attractive or repulse, 
weak or strong, short range or long range; all that matters is that it is Galilean invariant. One conse¬ 
quence of this invariance is Vij = Vji, so the two-body matrix elements also have the property 


V 


aj3 

"yS 


= V' 


^o 
5-y ’ 


(40) 


Putting these together, the following four two-body matrix elements are equivalent for any Galilean- 
invariant two-body interaction potential Vij: 


oB 


= V 


= V' 


^o 

5‘y 


= V 


6'y 


(41) 
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To better represent this symmetry, denote this two-body matrix element by vya-y\ yps\ i where the 
order or the symbols within the pair and the order of the pairs is arbitrary. This matrix element is 
proportional to the direct, first order transition amplitude for state transitions a <->■ 7 and /3 O 5. 
Additionally, for the contact interaction the two-body matrix elements also have the property 


olB cx. 3 Boi 

^5 = ^^57 = 


(42) 


because 


(a/3|Vi2|7(5) = g J dqidq 2 (j)a{qi)(l)p{q 2 )S{qi - q2)(l)j{qi)(t>5{q2) 

= 9 dq(j)a{q)4>fiiq)(l^'yiq)4’s{q)- 


(43) 


Combining (41) and (42), for contact interactions the two-body matrix elements with all of the 


24 possible permutations of the state labels take the same value. Therefore, for the contact interaction 
the notation vyap-ys\ is convenient. This state permutation invariance of the contact interaction is not 
shared by the Hamiltonians Hq or but it provides an alternate explanation for why the states in 
totally antisymmetric representations of Pat are unperturbed by the contact interaction. Since is 
symmetric under exchange of particles and the two-body matrix elements vy^p^gy are symmetric under 
exchange of states, annihilates states in the antisymmetric subspace K}d 1 . 


4.2 Conhguration Space Symmetry Group Cjv 

Any Galilean-invariant operator has the same configuration space symmetry as the coincidence 
manifold of all configurations in which at least two particles coincide in position 

N 

= (44) 

i<j 

where Vij is the (A^—l)-dimensional hyperplane with qi = qj. The coincidence manifold divides con¬ 
figuration space into N\ identical sector^ By Galilean invariance, the operator and manifold 
are invariant under permutation of particles and under total inversion 17. Additionally, Galilean 
® This is an essential difference between one dimension and higher dimensions: in higher dimensions the 
particles can slip past each other without the configuration passing through the coincidence manifold V^. 
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invariance implies that commutes with the total momentum P = ^ Pi, and therefore is invari¬ 
ant under translations in (or inversions of) the center-of-mass coordinate Q (x The configuration 

space symmetry of the coincidence manifold is therefore isomorphic to 

C]; ~ X 0(1) X (0(1) X Tr), (45) 

where are translations along the center-of-mass axis and the second copy of 0(1) is reflections 
Ur perpendicular to the center-of-mass axis. The other copy of 0(1) is total inversion 77. Relative 
inversion, i.e. inversion of the relative coordinates but preserving the orientation of the center-of-mass, 
is defined 77^ = UrII = PPIr and is also in the group C]^. 

Excluding the translation symmetry (which will be broken by any trap) the remaining symmetry 
group Sat x 0(1)^^ is isomorphic to the point group of a A^-dimensional prism with end faces that are 
dual 7V-simplices, e.g. a hexagonal prism for N = 3 (see previous article) or an octahedral prism for 
N = 4: (see Fig. [^. 

The configuration space symmetry group Cat of the total interacting Hamiltonian must contain 
the intersection of the transformations in C](- and C^. Consider the three cases: 

— Asymmetric well: Ca/ ~ Sat- The interacting Hamiltonian has the same configuration space sym¬ 
metry as the non-interacting Hamiltonian, namely just Sat. As described above, particle exchanges 
are realized in by orthogonal transformations that leave the center-of-mass axis and orientation 
invariant. 

— Symmetric well: Cat ~ Sat x 0(1). Total inversion 77 is a symmetry transformation and the 2N\ 
elements of Cat are realized by orthogonal transformations that leave the center-of-mass axis (but 
not necessarily the orientation) invariant. There is a doubling of irreps compared to the previous 
case, one set of [fj] G P{N) with even parity and one set with odd parity [/r]“. Note that even 
for symmetric traps relative parity is not a good quantum number for N > 2 interacting particles 
unless the external single-particle potential V^{Qi) is quadratic in position. 

— Harmonic well: Cat ~ Sat x 0(1)^^. For quadratic trapping potentials like the harmonic well, the 
center-of-mass and relative coordinates are separable. The total inversion 77 and center-of-mass 
inversion 11 r (or equivalently 11 r and relative inversion 77^) are independently good symmetry 
transformations. The group Sat x 0(1)^^ has order 4A^!. The four irreps [/r]^ for every partition 
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Fig. 1 This figure depicts the intersection of the four-dimensional, four-particle coincidence manifold (the 
six disk planes) with a three-dimensional hyperplane of constant center-of-mass gi -I- 52 + ?3 + 94 = c, i.e. 
the intersection of octahedral prism V"* and the relative configuration space. There are 24 sectors Qs for each 
s G S4 corresponding to a specific particle order 9 ^^ < qs^ < < Qs^- Reflections across the disk planes are 

the geometrical realizations of transpositions like (12). Three-cycles like (123) correspond to rotations by 2tt/3 
around the arrows (red online). Rotations by tt around the dashed line segments (blue online) correspond to 
two two-cycles like (12) (34). Rotoreflections by tt/2 around the solid line segments (green online) are four-cycles 
like (1234). 


[^] G P{N) correspond to the four possible combinations of total and relative parity. Note that 
for the harmonic well, compositions \v\ of non-interacting states have definite total parity but 
such compositions need not have definite relative parity (or center-of-mass quantum number, see 
below). There are multiple compositions with the same energy for the harmonic well, there are 
other observables that better interpolate between the U(iV) symmetry of the non-interacting case 
and the Sat x 0(1)^^ symmetry of the interacting case, namely the Jacobi hypercylindrical basis. 
For more details, see Ref. m- 
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4.3 Kinematic Symmetry Group Kjv 

Unless there are accidental or emergent symmetries, the kinematic symmetry group for interacting 
particles is just PArxKi~CArxTi, where in this case the time translation group generated by 
is just a phase in each energy eigenspace. The clocks of the individual particles are synchronized 
by the interaction and are no longer independent symmetries. For asymmetric traps and symmetric 
(non-harmonic) traps, the irrep structure of Kjv differs from Cat only by the addition of an extra 
quantum number, the total energy. 

For the harmonic trap, there is the additional kinematic symmetry factor U(l) due to the separa¬ 
bility of center-of-mass and relative coordinates. The center-of-mass degree of freedom in phase space 
is unperturbed by the interaction. The kinematic symmetry for the harmonic potential is therefore 
xO(l)xU(l)xTt. The U(l) symmetry gives each irrep of this group an additional quantum 
number n. There are an infinite number of irreps of Kjv labeled none of mix under interac¬ 

tions. For the case of contact interactions, this is exploited in Ref. [5^ [571 HH] to further decompose 
the kinematic Hilbert space into non-interacting subspaces that aid analysis and and efficient exact 
diagonalization and in Ref. |19j to construct adiabatic mappings. 


4.4 Approximation Methods 

Analytic solutions of the model Hamiltonian are only known for a few trap shapes and interactions: 

— harmonic traps with harmonic interaction^ 

— harmonic traps with contact interactions of any strength for two particles [551 EOl IHI] ; 

— hard wall traps with contact interactions of any strength for N particles (Bethe ansatz solution) [151 

naisiES]. 

For almost everything else, we need approximation methods. Symmetry makes it easier to make these 
approximations in an analytically and computationally efficient way, and provides spectroscopic clas¬ 
sifications that aid in the interpretation of the results. 

^ This system has kinematic symmetry Kat ~ U(l) x U(A — 1) x This is an example of an emergent 


interaction symmetry. 
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4-4-^ Exact diagonalization in a truncated Hilbert space 


As an example, consider exact diagonalization in the non-interacting basis. In Section 3, the energy 
eigenspaces /C of the non-interacting Hamiltonian were shown to carry irreps of and to be 
isomorphic to permutation modules The reduction of the the spaces into S^v irreps 

only depends on the shape of the composition [ly]. This critical observation is best exploited in the 
double tableau basis |ITH) for 1C. Because of the Sn symmetry, the matrix elements of the interaction 
are diagonal in Y: 

{WY\V^\W'Y') = {W\\V^ \\W')Syy', (46) 

where (W\\ \\W) is a reduced matrix element. This reduced matrix element is a linear combination 

of two-body matrix elements like ?^[a/ 3 j[ 7 ( 5 j for generic two-body interactions or for contact 

interactions. If the coefficients {n\WY) are known then one calculates 


(IT||y^||IH') = ^ (ITy|n)(n|H^|n')(n'|IT'y). 

nG .n'G [v^r\ 


(47) 


As an example, the reduced matrix element 

V2 


) IS 


(IP IT) ^ + n[^ 2 , 52 j - 2 n[^ 2 , 52 j) . 


The explicit expressions for (IT|| \\W') in terms of two-body matrix elements are sometimes unnec¬ 

essary for calculating physical results; see [iniiii] for examples of evaluating sum rules incorporating 
similar symmetrization techniques. 

Note that Y = Y' implies that all tableaux have the same shape [IT] = [IT'] = [T] = [T'], and so 
there are only matrix elements between vectors in the same tower subspace of irreps. This 

leads to two reductions in the computational scaling of exact diagonalization. First, one only needs 
to calculate matrix elements between states with the same shape [IT], and this means a reduction 
of necessary matrix elements by roughly a factor of the number of irreps for N. Second, for multi¬ 
dimensional irreps [/r], only one matrix element between basis vectors with the same Young tableaux 
Y is necessary (all the rest are the same), providing another reduction for each sector of (i[IT]. 

To see how this works in a truncated Hilbert space, consider four particles in an asymmetric trap 
and compositions \ i'\ such that 


0 • jzo + 1 ■ + 2 • 1^2 + 3 • 1/3 -I- 4 • z24 < 4 . 
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There are twelve compositions in this truncation: 

[O'*], [0^2], [O^aj, [0^12], [Ol^J, [0^4], LO^iaj, [0^2^], [01^2J, 

and there are 70 basis states in these twelve composition^ However, for spinless bosons, only the twelve 
symmetric states with [W] = [4] (one for each composition) are necessary for exact diagonalization in 
this truncated space. Spin-1/2 fermions could populate spatial states in subspaces K.^ with [W] = [2^], 
[21^], [1^], but there are no states with [W] = [1'*] in this truncated Hilbert space. The ground state for 

these fermions must have spin-0 and therefore only five states chosen from the [2^] subspace (different 

[X0\ 

W but all chosen with the same Y, e.g. ffl) are needed for exact diagonalization. 

The larger N, the more irreps there are, and for a given truncation dimension, the greater the basis 
reduction provided by symmetrization. However, the dimension of the truncation necessary to achieve 
a target accuracy grows more rapidly with N than the reduction due to symmetrization. So although 
symmetry methods reduce the computational challenge, the scaling of exact diagonalization remains a 
practical obstacle. 


4--4-2 Weak perturbations 


The same symmetry properties that facilitate exact diagonalization also assist degenerate perturbation 
theory for weak attractive or repulsive interactions. For example, unless there are emergent interaction 
symmetries, in first order perturbation theory each non-interacting subspace ^ splits into 


some number of Sat irreps [ji] determined by the Kostka number (c.f. eq. (lOl). For composition 

spaces that are simply reducible into Sn irreps, to first order in the weak perturbation the level 
with energy splits into levels carrying the irrep [/x] for each [/x] < \v\. The energy shifts depend on 
the trap shape and interaction, but the structure of level splitting and the first order energy eigenstates 
do not. For four particles, composition spaces like /CL""*!, and have this property. 

However, for composition spaces that are not simply reducible, the interaction breaks state per¬ 
mutation symmetry and the Weyl tableaux W no longer provide good quantum numbers for irreps 

that appear multiple times in the decomposition of K, . To find the first-order energy eigenvectors 
® Note that except for the harmonic trap, truncating the Hilbert space only to these compositions may or 


may not be consistent from a total energy perspective, but it is still an illustrative example. 






31 


using degenerate perturbation theory means diagonalizing the matrix formed by the reduced matrix 


elements (bb|| \\W) defined in (46 1 for all Vb, W in a composition space with the same shape. 

In the previous article, this done for iV = 3 for the mixed symmetry irrep [21] in the non-simply 
reducible compositions with shape \v\ = [1^]. For iV = 4, the compositions with shapes [21^] and [l^j 
are not simply reducible. For example, two copies of appear in K\°‘ .That means that in order 

to find the first order energy eigenvalues and eigenvectors of the two levels which carry the S 4 irrep 
[31] into which splits, a two-by-two matrix must be diagonalized: 

/, 


a 


a 


Using expressions for these matrix elements for the contact interaction, these two three-dimensional 
levels have the energies 


- + 5U|^q,2^2j -|- 2u|^^2^2j 


(48) 


± y^9(u|^Q,2^2j)2 - 14 u|^q,2^2jU|^q,2..j,2j + ^ {v --iv ^2 yV yp'i ^2 y - y -|-4(ti|^,g2.y2j )2 j 

This algebraic expression for the energy shift (and similar expressions for the first order energy eigen¬ 
states) manifests the weakest form of algebraic universality since the two-body matrix elements are 
explicitly required. For non-contact interactions, the matrix elements are more complicated and the 


algebraic expression is longer than (48), but has the same form. To solve level splitting in levels like 


jqVo‘M^\ ^ a three-by-three matrix must be diagonalized for irreps [31] and [21^]. That requires solving a 
cubic equation, so the algebraic expression are substantially longer, but still universal in this weakest 
sense. 

For N > 5, the Kostka number Ky^yy^y for a given composition can be greater than four, and 
therefore expressions for the first order for some irreps require solving a characteristic equations for 
which we have no a priori guarantee that a universal algebraic solution exists. While certainly there 
are algebraic solutions for some polynomials of order five and above (contact interactions in an infinite 
square well provides an example), I hypothesize that such algebraically-solvable high-order polynomials 
are idealized limiting cases and not typical for physically-realistic interacting systems. 

If there are additional conserved quantities, like parity for symmetric traps or center-of-mass exci¬ 
tation for harmonic traps, then the spatial Hilbert space JC can be further reduced into sectors that do 
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not mix under interactions. As before, this deeper reduction allows more numerically efficient numer¬ 
ical solutions. For the case of parity, there is a doubling of the number of towers one for each 

parity. However, there is no change to the results of first order perturbation theory because each com¬ 
position space has a parity determined by its composition. For harmonic traps, the separability 
of center-of-mass and relative degrees of freedom makes the double tableau basis less useful, because 
state permutation symmetry does not commute with the unitary group U(l). However, this additional 
complexity is more than compensated by the emergent kinematic symmetry Kjv SjvxO(l)xU(l)xTt. 
See [19] for some examples of how this scaling works for harmonic traps. Further techniques for max¬ 
imally exploiting this additional symmetry for harmonic traps are the subject of an article currently 
under preparation. 

5 Unitary Limit of Contact Interactions 

In the unitary limit g —>■ oo of the contact interaction in one-dimension, the particles cannot get past 
each other. Classically, the particles rattle back and forth in the trap, bouncing with perfectly elastic 
collisions off each other or rebounding from the edge of the trap potential. Quantum mechanically, the 
scattering from the delta-potential is diffractionless [33] and this system is called a Tonks-Girardeau 
gas. 

The order of the particles is stable under these dynamics. The coincidence manifold divides 
Q into N\ ordering sectors. Each order can be labels by a permutation s S Sat and corresponds to a 
sector of configuration space Qs C defined by the condition < qs^ < ■ ■ ■ < qs^- The sector Qs 
is bounded by the (N—l) hyperplanes Vsisai ..., and Vs„_isw See Fig. for a depiction of the 

four-particle coincidence manifold and sectors. 

In the unitary limit of the contact interaction, the wave functions for finite energy states must vanish 
at the edges of the sectors Qs, but inside the sector they satisfy the non-interacting Hamiltonian ■ 
Denote by L‘^{Qs) the Hilbert space of wave functions on Qs satisfying the nodal boundary conditions 
on V'^. All the ordering subspaces K-s ^ L^(Qs) are equivalent. Define the spatial Hilbert subspace 
/Coo by 

/Coo = 0 JCs. 


( 49 ) 
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What about states in /C but not in /Cqo? They have infinite energy because they have some probabil¬ 
ity density on or they have infinitely sharp discontinuities Note that states in a super-Tonks- 
Girardeau gas (corresponding to the unbound states in the limit g —oo) are also in /Coo, but the 
infinite tower of infinite negative-energy bounds states are not. 

The energy spectrum within each these independent sectors ICg can be deduced from the one-particle 
spectrum cti quite simply: whenever there is a fermionic state in the Wparticle non-interacting spec¬ 
trum there is a stationary state in JCg. This is because the boundary condition imposed by contact 
interactions is ‘automatically’ solved by the non-interacting fermionic states due to antisymmetry. At 
the unitary limit, wave functions must vanish on the coincidence manifold V'^. From a given one- 
particle spectrum ui, the spectrum and many features of the states at the unitary limit can be 
determined universally, including the incorporation of identical particles with spin 0 la [31HH [33 
13111371131133110 ]. All ordering sectors /Cg are identical and so the A^-particle spectrum cr“ in the 
unitary limit is composed of A^!-degenerate energy levels. 


5.1 Configuration Space Symmetries 

The configuration space symmetry of the Hamiltonian at the unitary limit must contain Cat 
as a subgroup because it is just a special case of a Galilean-invariant interaction. Therefore particle 
permutations p S Ptv are valid symmetries. The configuration space representation 0{p) maps sectors 
Qs onto each other like 

0{p)Qs = Qsp-i. (50) 

and the unitary representations map subspaces /Cg onto each other like 

17(p)/Cg = /Cgp-i. (51) 

These transformations are linear and continuous in . However, there are other configuration space 
transformations in that are non-linear and non-continuous, but whose unitary representations on 
® Technically, one can construct normalized superpositions of finite-energy states in /Coo that have expectation 
values for the total energy that diverge, but ignoring these kinds of pathological states 1 will refer to the space 
/Coo as the spatial Hilbert space of finite energy states. 
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/Coo commute with Namely, any other map that permutes the sector labels s like 

Ms = s' (52) 

is a configuration space symmetry. In this notation, particle permutations are represented by iV x TV 
matrices that act on the vector space of sector labels like 

M.{p)s = sp~^ 

and this is another example of the defining representation of Pat ~ Sat- 

The set of all maps like M form a group Qat that is isomorphic to the symmetric group S^v!- The 
particle permutations Pat C Qat are the only such maps that have continuous representations on Qat- 
All the other maps M G Qs are discontinuous, but because of the nodal structure V'^, the continuity 
of wave functions in /Coo is not disrupted. 

Ordering permutations Dn are a subgroup of sector permutations Qn that are complementary to 
Pat. Whereas p €Pn exchange particle labels (i.e. the numbers) in a sector s = {siS 2 • • • sat} wherever 
they occur, ordering permutations o G Dn exchange the order of labels in s, no matter what number 
is in that place in the order. Ordering permutations are indicated by letters, e.g. (ABC) means si 
becomes S 2 , S 2 becomes S 3 , and S 3 becomes si, or in other words they act like normal permutations 
on the sectors 

M(o)s = os. 

As an example with four particles, the particle permutation ( 12 ) exchanges the numbers 1 and 2 
in each sector s, e.g. M(12){3412} = {3421} and M(12){2431} = {1432}. In contrast, the ordering 
permutation (AB) exchanges the first and second number in the permutation s, no matter what the 
numbers are, e.g. M(Ai?){3412} = {4312} and M(Ai?){2431} = {4231}. Therefore, in addition to 
carrying a representation of the normal particle permutation symmetry, /Coo also carries a copy of 
ordering permutation symmetry D n that (like symbol permutation symmetry in permutation modules 
is isomorphic to Sat- One way to denote the configuration space symmetry group of is 
therefore 

D Piv X Dat X Cl - ^ Cl. 

A similar construction holds for 


(53) 
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The rest of this section applies these dual symmetries of particle permutation and ordering per¬ 
mutation. The inclusion of parity and additional kinematic and dynamic symmetries is also discussed. 
The combination of these symmetries provide an alternate, complementary method for treating the 
near-unitary limits as those described in 


5.2 Snippet Basis 


qlMr];s) = <^ 


One basis for JCa^ is provided by the snippet basis [TJ IMl m]. Denote by | [i/J [1^]) the single, totally 
antisymmetric state in a composition space that has composition shape [v\ = [1'^]. Then the 
snippet basis vectors are the A^! states denoted | \v\ s) with the property 

(54) 

0 qS Qs 

Each spatial Hilbert space sector K-s is spanned by the infinite tower of states for all 

[v] = [1^]. The A^! snippet basis vectors span an A^!-degenerate energy level with 

energy . 

Sector permutations M G Q 3 transform snippet vectors as expected 


and specifically particle permutations p G Ptv and ordering permutations 0 G Dat have the represen¬ 
tations 


17(p)|H[l^];s) = |H[l^];sp-i)and 
17(o)|H[l^];s) = |H[l^];os). 


(55) 


From this representation it is clear that particle permutations and ordering permutations commute. 


The representations (55) of both copies of Sat are A^!-dimensional and therefore must be reducible. 


In fact they are both the regular representations of Sat and the representation space is isomorphic 
to the permutation module 1. As with state permutation symmetry in compositions spaces 
with [v] = [ 1 ^], the canonical subgroup chains Sat A Sat-i D • • • D S2 and Dat 3 ‘^n-i D • ■ ■ D D2 
can be used to construct a complete set of commuting observables for the space . The eigenvectors 
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for these operators can be chosen as a double tableau basis F), where is a Young tableau 

filled with symbols A, B, C, etc. that denote the irrep in which the ordering permutation subgroup 
chain is diagonalized and the Young tableau Y filled with symbols 1, 2, 3, etc. denotes the irrep in 
which the particle permutation subgroup chain is diagonalized (see also Ref. 0 )- This notation was 
used for fV = 3 in the previous article. For iV = 4, the twenty-four basis vectors constructed from the 
fermionic state in are 


fijSgESEEES ) = 
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In this basis, the S 4 irreps given by the shape [S^] = [F]. The coefficients between the snippet basis 
and the double tableau basis ([z^J [1'^]; s| [1-^]; F) can be found efficiently using computer-assisted 
algebra and the methods of Ref. |3] . 

Note that every composition \ v\ with arbitrary shape [ly] is associated with another composition 
[iy'\ with shape [1-^] by adding the sequence (0,1, 2,..., (TV—1)) to the representative ordered sequence 
n € [z^J. Therefore there is also an element | \ v'\ s) G JCp for every composition [z^J, i.e. the Bose- 
Fermi mapping ED- 


5.3 Incorporating Additional Trap Symmetries 

One element Ojr G Oat of the ordering permutation group is equivalent to reversing the order of the 
particles. For example, for three particles o,r is (AC) and for four particles Ojr = {AD)(BC). For 
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symmetric wells, the parity operator is realized on the snippet basis using this operator: 

n\ H [1^]; s) = n\ H [l^]; {siS 2 • • • sn-iSn}) 

= ILH [1^]; {snsn -1 ■ ■ ■ S 2 S 1 }), (56) 

where ’ ’ ’ '^I'r is the total parity of the composition [z/J. This representation of the parity 

operator can be diagonalized to decompose into Sw x 0(1) irreps . These will not necessarily 

be the same states constructed using the order permutation symmetry subgroup chain that leads to the 
double tableau basis because the operator 0 ^ is not part of the subgroup chain Dn 3 £ 1 az-i D • • • 02 - 
However, 0 ^ can still be diagonalized along with the normal particle permutation group subgroup 
chain, and an alternate ordering permutation subgroup chain can be found. For example, one can show 
that for four particles a complete set of commuting observables that respects parity is 

^U234J ^ ^^^ 2 ) + f/(i3) + (7(23) + 17(14) + 17(24) + 17(34), = 17(12) + 17(13) + 17(23), 

^ y(^ 2 ), 17(o^) = U{{AD){BC)), and = U{AD). (57) 

This notation for conjugacy class operators is a slight modification of Ref. [3], and it is extensible to 
larger N. 

Not all of the Pat ~ Sat irreps have the same parity as the parity of the original fermionic com¬ 
position. For example, the subspace with = 1 reduces into the following S 4 x 0(1) 

irreps 

/C ^ ^[4]+ _^[ 3 i]+ 33 0 33 33 2 ^ 1311 - 33 2^121^]-. (58) 

Looking at the dimensions d[fj], we see that half the states have positive parity and half have negative 
parity; this result holds for all N. 

Besides the examples with N = 2 and N = 3 in the previous article, further examples of this 

diagonalization are discussed for the harmonic case in [T^] and tables for iV = 3, N = 4 and N = 5 

are provided. In the harmonic case, relative parity and the center-of-mass excitation are also good 
quantum numbers. They are not necessarily commensurate with composition subspaces, but they can 
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be diagonalized simultaneously with the particle and ordering permutations symmetries. The extra 
symmetry means that the unitary limit in a harmonic trap is superintegrable, and therefore in the 
near unitary limit I conjecture the system is likely to resist thermalization longer than less symmetric 
traps. 


5.4 Near Unitary Limit 

In the near unitary limit, ordering permutation symmetry is broken because there is tunneling between 
adjacent sectors Qs- At unitarity, the kinematic symmetry of the Hamiltonian is D Pn x x 
K i, but away from unitarity the symmetry is only K^v Pw x Ki. Therefore, each energy level 
splits into the Pat irreps that compose it. 

The tunneling can be represented by an operator that breaks the ordering permutation symmetry 
and ‘lowers the energy cost’ of having a cusp on the coincidence manifold V^. Then near unitary the 
Hamiltonian is split perturbatively into + T. A suitable operator T has the form 

T = —aAsU^AB) — ascU{BC) — ■ ■ ■ — aMNU{MN) — {gab + asc H— • )U{e), (59) 


where the real, positive coefficients gab, o,bc^ etc. are the tunneling amplitudes for adjacent sectors m- 




= — [ d^(iS{qi-qj) 
9 Jo., 


a(q|Hr]) 


dqt 


(60) 


The last term in (591 that is proportional to the identity renormalizes the energy shift so that the 


totally-antisymmetric spatial state feels no effect of the tunneling. If the trap is symmetric, then 
o-AB = clmNi olbc = o-LM, ctc. For understanding level splitting and the formation of bands that 
depend on trap shape, the ratio of the tunneling amplitudes (and not the absolute scale) is important. 

For example, see Fig. (also reproduced in the introduction of the previous article). With four 
particles in a symmetric trap, the tunneling operator has the form 


f = -tU{AB) - uU(BC) - tU(CD) - {2t + u)U{e). 


(61) 
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Fig. 2 Approximate level splitting diagram for four particles in three symmetric traps: (a) double well; (b) 
infinite square well; (c) V-shaped or cusped well, i.e. softer than harmonic. The thick band on the right in each 
figure is the 24-fold degenerate ground state energy level for four distinguishable particles in the unitary limit 
of the contact interaction. The three trap shapes are distinguished by the ratio of tunneling amplitudes tju. 
The following ratios have been chosen to illustrate the trap dependence of these amplitudes: (a) tju = 2.9; 
(b) tju = 1; (c) tju = 0.3. The idea is that (a) for double wells tunneling in the middle is suppressed so 
t > u; (b) for inhnite square well the potential is uniform inside the trap and so (for low particle density) t 
and u are approximately the same; (c) for softer wells, there is more phase space in the middle of the well so 
u> t. For harmonic wells, t/u ~ 0.762 (c.f. [411142114^ I44j l. In subfigure (c) a more extreme ratio is depicted, 
corresponding to a V-shaped or cusped trap. 


Diagonalizing the basis provided by (57) shows that the energy splits into an energy level with t^Yv\ =1 


for each S 4 x 0(1) irrep in the sum (58) as follows: 


[4]+: 

—At — 2u 

[31]+ : 

-2t - 2u 

[ 21 + : 

-2t-u± 

[ 2 l 1 + : 

-2t 
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—3f — M ± \/t^ - 

[ 2 ir : 

—t — + 


( 62 ) 
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Unlike the three particle case, parity does not completely remove the trap dependence of the first 
order energy eigenstates; there are two copies of [2^] + , [31]“, and [21^]“ that require an additional, 
trap dependent operator to diagonalize. 

If instead an asymmetric trap was considered, there are three tunneling parameters for U{AB), 
U{BC) and U{CD). There are three copies of irreps [31] and [21^], so computing level splitting requires 
solving a cubic equation. For = 5 in asymmetric traps or iV = 6 in symmetric traps, the required 
diagonalizations require solving greater than quartic equations, and the last traces of algebraic uni¬ 
versality are lost for the near-unitary regime. The characteristic polynomials have algebraic solutions 
for certain trap shapes (like the infinite square well), but as with the case of the weak level splitting, 
I hypothesize that those cases are not typical. 

When the level spitting problem is algebraically solvable in the weak interaction limit and the near 
unitary limit, then the adiabatic mapping problem between non-interacting states and their unitary 
limit is in principle solved. The order of the levels within a particle tower should not change unless 
for some particular value of (/ a new multiparticle symmetry emerges. That does not seem likely for 
the contact interaction, but I do not have a proof that it is impossible. 


6 Conclusion 

The essential claim of this pair of articles is that the configuration space and kinematic symmetries 
of Hq , , and provide a powerful and unifying tool for analysis and computation. One way to 

specify the focus of this article is to ask the series of questions: 

— Given only the symmetries, what can we say about the A^-particle spectrum apfl How much about 
the energies, degeneracies and states for the interacting system can be inferred without any specific 
knowledge of the one-particle spectrum a±7 How does the spectrum change when the interaction 
strength is tuned adiabatically or rapidly quenched? 

— If we also know the specific spectrum cti in addition to the symmetries, how much more can we say 
about the spectrum ctat? And how much more if we also know the wave functions ipniQ) = {qW) of 
the one-particle energy eigenstates H^\n) = e„|n)? 
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— Finally, what if we know the explicit form of the two-body interaction Vij and/or can calculate the 
two-body matrix elements {nin 2 \Vi 2 \n'in' 2 ) = L"'i" 2 J between noninteracting states? 


As a general conceptual framework, the less we have to know about the specifics of the trap or interac¬ 
tion, the more ‘universal’ the results. What this article demonstrates is how this notion of universality 
breaks down for most interacting system as the number of particles is increased and the symmetries 
of the trap and interaction are reduced. 

Despite its length, this article has left out many relevant topics, like the SO(2,1) dynamical symme¬ 
try of the contact potential in a harmonic trap, lattice symmetries, supersymmetric potentials in one 
dimensions, and the possibility of interaction symmetries that depend on the internal structure of the 
particles. Further, these symmetry classihcations could be generalized to higher dimensions. Symmetry 
should be less constraining as the number of degrees of freedom grows and algebraic universality may 
break down sooner. Also, the effect of intrinsic three or higher few-body interactions on the spectrum 
would be interesting to incorporate into symmetry analysis to seek possible manifestation in spectral 
shifts. 

However, without adding to the complication of the one-dimensional trap, two-body interaction 
model, there is still much work to be done. Efficient methods of state construction are required for 
perturbation theory and exact diagonalization and for calculating reduced density matrices, correlation 
functions, and entanglement spectra among particles and between spin and spatial observables. Another 
possible avenue for future work is to exploit the close connection between finite groups and number 
theory. Perhaps there are practical protocols for simulations of number theory problems that employ 
combinations of adiabatic tunings and diabatic quenches of the trap shape and and interaction strength 
to manipulate states. 
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Table 1 This table provides information about the kinematic symmetry group and configuration space 
symmetry group Cjv for N = 2, S and 4 identical non-interacting particles in asymmetric, symmetry and 
harmonic traps. See text to understand notation for structure of Cjv and K^v- For each group, the dimensions 
of the unitary irreducible representations (irreps) is provided, and the number of inequivalent irreps with that 
dimension is indicated by the superscript. The symbol A £ N labels inequivalent 0{N) irreps and X £ N 
labels inequivalent U(A'^) irreps. Note that for N = 2, there is an infinite tower of inequivalent irreps of 0(2) 
for A > 0. The Coxeter notation and order is given for the finite-order configuration space groups, and the 
Schonfiies point group notation is also given for N — 2 and 3. 



Asym. Trap 

Sym. Trap 

Harm. Trap 

o 

too 

Sa 

SaxO(l)^'' 

0(2) 

Irreps 


1^,2 

1,2 

Point 

Di 

D 4 

0(2) 

Coxeter 

Ai~[] 

BCa ~ [4] 

— 

Order 

2 

8 

CX) 



Sax(0(l)xTt)^2 

U(2) 

Irreps 

1,2 

l2 23 

(^ + 1) 

c§ 

Ss 

S3xO(1)><3 

0(3) 

Irreps 

1^2 

1‘‘,2^3‘‘ 

2A-f 1 

Point 

C3„ 

Oh 

0(3) 

Coxeter 

Aa ~ [3] 

BC 3 ~ [4, 3] 

— 

Order 

6 

48 

CX) 

K§ 


S3x(0(l)xTt)^3 

U(3) 

Irreps 

1,3,6 

1^3^6‘‘ 

(X+2Y. 

2X\ 

Cl 

S4 

84x0(1)'' 

0(4) 

Irreps 

1^2,3^ 

1^2^3^4^6^82 

(A + 1)2 

Coxeter 

A 3 ~ [3, 3] 

BC 4 ~ [4, 3, 3] 

— 

Order 

24 

384 

CX 

K5 

S 4 xTt’'‘‘ 

S4x(0(l)xTt)>'4 

U(4) 

Irreps 

1,4,6,12,24 

1^4^^, 6®, 12®, 24® 

(X-l-3)! 

6X! 
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Table 2 This table provides information about the kinematic symmetry group Kjv and configuration space 
symmetry group Cjv for N = 2, 3 and 4 identical non-interacting particles in asymmetric, symmetry and 
harmonic traps. For each group, the dimensions of the unitary irreducible representations (irreps) is provided, 
and the number of inequivalent irreps with that dimension is indicated by the superscript. The Coxeter notation 
and order is given for the finite configuration space groups, and the point group notation is also given for N = 2 
and 3. 



Asym. Trap 

Sym. Trap 

Harm. Trap 

Ca 

Sa 

SaxO(l) 

SaxO(l) 

Irreps 




Point 

Di 

Da 

Da 

Coxeter 

Ai~[] 

Ai^ ~ [2] 

Ai^ ~ [2] 

Order 

2 

4 

4 

Ka 

SaxTt 

SaxO(l)x Tt 

SaxU(l)xTt 

Irreps 




Ca 

Sa 

SaxO(l) 

SaxO(l)^2 

Irreps 

1^2 

1\2^ 

CO 

Point 

Ca„ 

Dad 

Deh 

Coxeter 

Aa ~ [3] 

[[3]] 

[[3], 2] 

Order 

6 

12 

24 

Ks 

SaxTt 

SaxO(l)x Tt 

S3xO(l)xU(l)xTt 

Irreps 

1^2 

1\2^ 

18,24 

C4 

S4 

S4XO(l) 

S4xO(1)^2 

Irreps 

1^2,3^ 

1^,2^3^ 

1®, 2'*, 3® 

Coxeter 

Aa ~ [3, 3] 

[[3,3]] 

[[3, 3], 2] 

Order 

24 

48 

96 

K4 

S4xTt 

S4xO(l)xTt 

S4xO(l)xU(l)xTt 

Irreps 

1^2,3^ 

1^,2^3^ 

1®, 2'*, 3® 




